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Abstract — Code annealing, a new method of designing good 
codes of short block length, is proposed, which is then con- 
catenated with cyclic lifting to create finite codes of low frame 
error rate (FER) error floors without performance outliers. The 
stopping set analysis is performed on the cyclically lifted code 
ensemble assuming uniformly random lifting sequences, and the 
suppressing effect/weight of the cyclic lifting is identified for the 
first time, based on which the ensemble FER error floor can be 
analytically determined and a scaling law is derived. Both the 
first-order and high-order suppressing effects are discussed and 
quantified by different methods including the explicit expression, 
an algorithmic upper bound, and an algebraic lower bound. 

The mismatch between the suppressing weight and the stop- 
ping distances explains the dramatic performance discrepancy 
among different cyclically lifted codes when the underlying base 
codes have degree 2 variable nodes or not. For the former case, 
a degree augmentation method is further introduced to mitigate 
this metric mismatch, and a systematic method of constructing 
irregular codes of low FER error floors is presented. Both regular 
and irregular codes of very low FER error floors are reported, 
for which the improvement factor ranges from 10 6 -10 4 when 
compared to the classic graph-based code ensembles. 

I. Introduction & The Main Results 

The asymptotic analysis of low-density parity-check 
(LDPC) codes perfectly describes the performance when code- 
word length n > 10 7 . Nonetheless, practical implementation 
generally requires short-to-moderate n < 10 4 and efficient 
code representation. Two major approaches of designing good 
finite codes consist of systematically searching through the 
random code ensemble or directly constructing codes under 
combinatorial/algebraic guidelines [1], [2]. The former ap- 
proach is best represented by the progressive edge growth 
(PEG) construction, in which greedy/trellis-based algorithms 
for edge assignment have been proposed with different ob- 
jective functions such as the girth [3], the extrinsic cycle 
degree [4], partial elimination of stopping sets [5], further 
criteria regarding degree 2 variable nodes [6], and the averaged 
error rate upper bound [7]. 

If one focuses on a finite code ensemble instead of a 
single code, recent work finds progress on different subjects 
including the error rate curves [8], the scaling law of the 
waterfall region [9], and the error- floor upper bounds [10]. 
As most researches along this direction, our work is based 
on the binary erasure channel (BEC). Existing simulations 
have demonstrated high correlation between the performance 
on BECs and on other channel models. 



This paper first discusses a novel code annealing (CA) 
algorithm for finite code optimization based on the tool 
developed in [11] and then studies the interaction between 
the base code and its corresponding cyclically lifted (CL) 
code ensembles. A new connection between the finite code 
optimization and the finite-sized code ensemble analysis is 
established by considering the "suppressing effect" of cyclic 
lifting. The CA algorithm uses any existing code as a starting 
point and gradually enhances its performance, and thus can 
be seamlessly combined with any existing methods. CA is 
also the first efficient code optimization with the objective 
function being the stopping distance (D stp ), defined as the 
minimum size of stopping sets, instead of approximations or 
design heuristics. No constraint on the number of degree 2 
variable nodes is assumed for CA and very good frame error 
rate (FER) error floor is reported for n = 576. 

For the CL ensembles, we show that the probability that a 
stopping set in the base code survives after lifting is 

where K is the lifting factor, #U is the size of the base 
stopping set, #E is the number of edges in the subgraph 
induced by the base stopping set, and #C dd is the number 
of corresponding check nodes with odd degrees. This formula 
prompts a definition of the suppressing weight W sup = 
0.5#i? — jfV + 0.5#C o dd, and designing a good base code 
(suitable for cyclic lifting) is equivalent to maximizing the 
Wgup of all small stopping sets of the base code. CA is 
capable of constructing a n = 128 regular (3,6) codes of 
minWgup > 5 so that all its stopping sets are experiencing 
at least K~ 5 suppressing effects. For K = 32, the lifted 
codeword length being n = 4096, the suppressing factor is 
2.98 x 10~ 8 and results in FER al)g w 7 x 10~ 14 at erasure 
probability e = 0.3. That is, with very high probability, a 
typical code will have extremely low FER, which demonstrates 
a perfect match between CA and the CL ensemble. With base 
codes being optimized by CA, the explicitly constructed CL 
ensemble has similar behavior as the expurgated ensemble 
[10]. 

The suppressing effect analysis also sheds light on designing 
irregular codes of many degree 2 variable nodes, which 
generally has min JV sup = or 1 and does not generate 
satisfactory results even when codes of large stopping dis- 
tances -D s tp are used as base codes. A combination of CA and 



degree augmentation (DA) is introduced to mitigate this metric 
mismatch between W sup and D stp , which is the first definite 
method of "neighborhood optimization" [12]. An n = 576, 
rate 1/2 irregular code with 62.5% degree 2 variable nodes is 
constructed with intact threshold performance and error floor 
FER w 4 x 10~ 5 at e = 0.3. This result shows that good error 
floors can still be obtained even in the case that there are many 
cycles involving only degree 2 variable nodes. 

II. Code Annealing (CA) 

A. Stopping Sets & the Asymptote of Error Floors 

For a LDPC code of size n and mxn parity checkmatrix H, 
we denote the coded symbols and the stopping 

set as s C {x±, ■ • • , x n } [8]. The error floor is determined by 
the minimum stopping sets, the size of which is defined as the 
stopping distance D stp , the order of the error floor, while the 
total number of which is termed as M s , the multiplicity of the 
error floor. M s x e D=lp then becomes an asymptote of the error 
floor. In [11], an efficient algorithm is provided to determine 
D s tp, Ms, and {s : |s| = D stp } for codes of n < 500. 

For a code ensemble C = {C a } a where C a represents 
the element code indexed by a, we define the order and 
multiplicity of the error floor of C by 

•Dstp.c = mm{D stp (C a ) : VC Q e C} 
M sfi = E{|{s(C a ):|s(C a )|=£) s tp, c }|}, (1) 

where D stp (C a ) and s(C Q ) are the stopping distance and 
stopping sets of C a . The asymptote of the ensemble FER 
error floor becomes M S) c x e Datp ' c . One example is the classic 
regular (3,6) code ensemble of n — 512, for which we have 

Atp.C = 2 and M sfi a 6.8 x lO" 2 . 

B. The Algorithm 

The CA algorithm is motivated by the following observa- 
tion, and similar phenomena were observed in [13]. 

• For arbitrary fixed codes, the FER in the error floor region 
is generally dominated by only a few "bad bits," while 
other bits are experiencing almost error free transmission. 
This typical uneven protection prompts room for improvement 
when properly utilized. By locally changing the intercon- 
nection of those "bad bits," the corresponding error patterns 
can be removed and the error floor performance can be 
drastically improved. A detailed description of CA is provided 
in Algorithm 1, in which the stopping set exhaustion method 
in [11] is adopted. This new technique is termed as code 
annealing, since most of the original code structure is kept 
unchanged while a small portion of edges is rearranged to 
generate a code of much higher performance. 

The advantages of CA are three-fold. First, only local inter- 
connection is changed during CA, and the waterfall threshold 
thus remains identical to its original value. Second, the local 
improvement nature makes CA compatible to all existing good 
codes or code ensembles, and can be used to construct codes 
from scratch or to polish existing codes. Another advantage is 
the use of the sufficient metric (Astp, M s ) instead of indirect 



Algorithm 1 The Code Annealing (CA) Algorithm 

1: Set the target order d = 1 
2: repeat 

3: Construct S, a collection of all stopping sets of size d. 
4: while S is not empty do 

5: Sequentially select s, x a , and e a as follows. Uniformly randomly 
select s from S, x a from s, and e a = (x a ,y a ) from all edges 
connecting x. 

6: Uniformly randomly select xt, from var. nodes other than s, select 

e b — ( x b>yb) from all edges connecting Xf,. 
7: Remove e a and ej,, and add e' a = (x a ,yb) and ej, = (zb,Ua) to 

the original bipartite graph. 
8: Find D s tp and M s of the new graph. 
9: if (D stp ,A/ s ) >~t (d, |S|) then 

10: Update S so that it contains all stopping sets of the new graph of 

size d. 
1 1 : else 

12: Remove e' a and e' b , and add back e a and e^,. 

13: end if 
14: end while 

15: d^d+1 

16: until the preset computation time limit is reached. 

f (di, mi) >- ((^2, m2) iff either d± > efe, or d\ — d2 and mi < 7712. 



ones such as the girth, the extrinsic cycle degree, etc. With this 
matched objective function, the benefit of CA is guaranteed 
and the amount of improvement is contingent only upon the 
efficiency of the stopping set exhaustion method. Namely, even 
with the random selection involved in the CA construction, it 
is guaranteed that there is no performance outlier in the code 
ensemble constructed by CA. 

The edge swap trials in CA also admit parallel computing, 
which makes CA an appealing choice compared to other 
algorithms for which sequential computation is necessary. 

C. Performance Comparison — Three Examples 

Example 1: A fixed regular (3,6) code of size 64, labelled 
as C\, is randomly chosen from the graph-based ensemble 
and has D stp = 4 and AI S = 1. After rearranging 41 out of 
the total 192 edges, the resulting code C*2 has D stp = 8 and 
M s = 89. The bit-wise error floor order is now 8 for all bits 
and even error protection is achieved. 

Note: For any regular (3,6) code of size 64, the maximum 
girth is < 3 (counting only the variable nodes) while D stp 
can be improved to at least 8, as demonstrated in C2. This 
highlights the importance of using the sufficient metric D stp as 
the objective function especially for short code optimization. 

Example 2: C3, a fixed, rate 1/2, irregular code of size 72, is 
randomly chosen from the graph-based ensemble with degree 
distributions [8], \{x) = 0.4187a; + 0.1626a; 2 + 0.4187x 5 and 
p(x) = x 5 , optimized for BECs. C3 has D s t p = 2 and il/ s = 
1. After rearranging 61 out of the total 216 edges, the new 
code C4 has D stp = 8, M s = 71. In C4, all bits except one 
have error floor order 8 and the exception bit has order 9. The 
error protection is again balanced. 

Example 3: A rate 1/2 irregular code of n = 576 is 
drawn from the same degree distributions as in Example 2. 
(D stp , M s ) is improved from (2, 2) to (11, 58) after rearrang- 
ing 247 out of the total 1726 edges. A significantly lower error 
floor is observed in Fig. 1(c) and the threshold performance is 




intact, which is best seen by the corresponding bit error rate 
(BER) curves. 

As plotted in Fig. l(a-c), the error floors of both the FER 
and the BER curves are significantly lowered by CA. In sum, 
the local edge optimization of CA is able to locate a very good 
code in the neighborhood of the starting point. 

III. The Suppressing Effect of the CL Ensemble 

A. Preliminary — Cyclic Lifting 

The idea of graph lifting can be traced back to [14] and 
is recently utilized in many coding researches [12]. The 
cyclically lifted (CL) code ensemble can be constructed by 
a base parity check code (or equivalently by a base parity 
check matrix H), a lifting factor K 6 N, and randomized 
lifting sequences 1. To be more explicit, we first index the 
l's in H by li, I2, • • • , where n e is the total number 
of l's. Uniformly randomly select a lifting sequence 1 = 

(h,h,-'- ,l ne ) from the set {°i !> ••• > K ~ 1 } n "- Let r ' 
denote a K x K permutation matrix obtained from circularly 

shifting an identity matrix I a distance I to the left. The lifted 

code based on the lifting sequence 1 then has a parity check 

matrix H' of size Km x Kn, which is obtained by replacing 

all O's by a K x K zero matrix and by replacing 1, by T li for 

all i S {1, • • • , n e }. It has been shown that the lifted structure 

admits efficient representation and encoding methods. 

B. The First Order Suppressing Effect 

Let Cl denote the lifted code obtained from a base code 
Cb- For any stopping set of the lifted code Cl, rename all 
Xi e Si by sGi<-ij +1 and use s L / K to denote the result after 
this conversion. "We construct s b from s L / K by removing the 
repeated nodes in s L / K . We then have 

Lemma 1: sb is also a stopping set of the base code Cb. 
In words, every stopping set in the lifted code corresponds 
to a unique stopping set in the base code. We then have the 
following corollary: 

Corollary 1: The sopping distance of the lifted code is no 
less than that of the base code. 

Based on Lemma 7, instead of selecting sl first as described 
in the previous construction, we sometimes directly choose sl? 
as a stopping set of Cb and define sj, as a stopping set of 



Cl that is compatible with s^. The suppression of different 
orders is then defined as follows. 

Definition 1 (Suppression of Different Orders): If $l/k 
contains no repeated nodes, we define sl as the survival of 
sb against the first order suppression. Sometimes we use "s l 
is a first order survival of sl?" as shorthand. The high order 
survivals are defined as those sl not a first order survival. 

With a fixed lifting factor K and a uniformly randomly 

chosen lifting sequence, we then quantify the first order 

suppressing effect in the following theorem. 

Theorem 1: For any s^, let Sl,i denote the collection of 
all first order survivals of s^. We then have 



E{|Si,i|} 

K #v-#E g( ' 3 "(-!)' ( deg t (Cj) ) ( K t )t\(K-t) d ^ 



#C ( mm[K, deg(o 3 )) 
j=l \ t=0 

where #V = \sb\, #£' and #C are the numbers of check 
nodes and edges in the subgraph (of Cb) induced by Sb, and 
Cj and deg(cj) represent the individual check node and its 
degree in the induced subgraph. For future reference, we define 
E{|Sl i|} as the first order suppressing effect of s^?. 
Proposition 1: For any s^?, E{|Sl,i|} is of the order 

^-(0.5#-E-#V+0.5#C od d) 

where the new notation #C dd is the number of check nodes 
of odd degrees in the subgraph induced by s^. 

Theorem 1 and Proposition 1 are derived from the enumer- 
ator polynomials for the CL ensemble. 

Define the suppressing weight of sb as W^up := 0.5#i? — 
#V + 0.5#C o dd- From the above results, the larger Ws Up is, 
the smaller expected number of the first order survivals will 
be after cyclic lifting. Furthermore, when combined with (1), 
one can obtain the following error floor scaling law. 

Theorem 2: For the cyclically lifted ensemble Ccl, we have 

Artp.CcL = D stp (C B ) 

M SjCcl = K-n^Wnp: all min. s B } . Ms (C B ){l + o(l)). 
For base codes of minimum variable node degree > 3, W sup > 
0.5D stp for all s B . Using a n = 128, (D stp ,M s ) = (10,40), 
regular (3,6) code obtained from CA, the ensemble averaged 
FER error floor can be significantly lower even for small 
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(a) The Graph Ensemble (b)The CL Ensemble 

Fig. 2. Comparison between the classic graph-based ensemble and the CL 
ensemble with K = 4, base code (-Dstp, M B ) = (10,40). The codeword 
length for both ensembles is n = 512, and 10 codes are randomly chosen 
from each ensemble. 

lifting factors K = 4-8 (Kn = 512-1024). The comparison 
of ten n = 512 sample codes chosen from the CL and the 
classic ensembles is provided in Fig. 2. With a much lower 
averaged FER error floor, (D stp £ CIj , Ms.Ccl) = (10,0.015), 
situated at 8.74 x 10~ 8 when e = 0.3, the performance of 
the CL ensemble is highly concentrated when FER > 10 -6 , 
and no error floor can be observed by the Monte-Carlo 
(MC) simulation. Without those apparent performance outliers 
as in the classic ensemble, the CL ensemble behaves very 
much like the expurgated ensemble. For comparison, under 
MC simulations, all 10 randomly sampled CL codes are of 
performance indistinguishable from the best existing n = 504 
regular codes optimized by PEG [3]. 

From our experiments, almost all of sizes close to 
£ ) stp(Ci) are first order survivals, which can be explained by 
the following two reasons. By definition, high-order survivals 
have repeated nodes in s L / K . Therefore, for any the 
sizes of high order survivals are strictly larger than the size 
of the first order survivals if there is any, and -D s tp(Ci) is 
more likely to be determined by the first order survivals. 
Secondly, high order suppressing effects are decided not only 
by combinatorial arguments as in Theorem 1 but also by 
algebraic constraints as in Theorems 3 and 4, which further 
limits the chances of seeing high-order survivals of small sizes. 

C. The High Order Suppressing Effect 

The expression to quantize the high order suppressing 
effect is generally more involved. Let #Cl, and #£l 

represent the number of variable nodes, check nodes, and 
edges in the subgraph of Cl induced by s^, and define #Vb, 
#Cb, and similarly. We then have 

Theorem 3 (Algebraic Lower Bound): For any s^, let 
Sl ; >i be the collection of sharing the same #Vt, #Cl, 
and #i?L values. It is either Sl ; >i is empty or 

E{|5 L ,>i|} > 

const • ma ^ K -{#E L -#v L -#c L )^ K -(#E B -#v B -#c B)) 



for some const not depending on K. 

Sl,>i contains either first order survivals or high order 
survivals exclusively, depending on whether = #Vb 

or not. Theorem 3 is derived from the linear homogeneous 
equations associated with stopping set configurations such 
that the high order survivals correspond to the number of 
non-zero solutions. In many occasions when the number of 
equations is larger than the number of free variables, Sl ; >i 
is empty, which significantly limits the high order suppressing 
effect. Theorem 3 is tight when considering the first order 
suppressing effect for which #Vl = #Vb, #£x = #Eb, 
and #C L < 0.5#£ B - 0.5#C7 B , odd . 

For any s^, the expected number of high order survivals 
can be determined by directly solving the rank of the linear 
equations imposed by s^. This computation however has to be 
repeated when different sg are considered. A simpler upper 
bound will be provided, which is applicable to all possible 
choices of sb- 

Let the repetition R(xi) of any xi G sb denote the number 
of repeated appearances of Xi in Sl/k- Let x .d. C sb denote 
a subset of variable nodes such that by revealing the values 
of all elements in x o c i., the rest variable nodes sb\xva. 
can be successfully decoded in a non-increasing order of the 
repetition. Namely, if xi is the next to be decoded variable 
node via check node y,j, there exists at least another variable 
node x^ ^ x.i connected to yj such that R(x' i ) > R(xi). The 
subscript "o.d." stands for ordered decoding. We then have the 
following theorem. 

Theorem 4 (An Algorithmic Upper Bound): If we slightly 
abuse the notation and use Sl,>i to denote the collection of 
sl compatible to a given repetition pattern {R(x)} x ^ Sb , there 
exists a const not depending on K such that Vif, Vx .d., 

E{|S L .>!|} < const- K^—^ {R{x) ~ 1] )k- w ^. 
Similarly, Sl.>i contains either first order survivals or high 
order survivals exclusively, depending on whether {R(x)} are 
all l's. Theorem 4 is tight both for the first order survivals 
and for the parallel high order survivals with |x o c i.| = 1 and 
R(x) — R,Vx £ sg. Theorem 4 further confirms the role of 
W sup for high-order suppression. 

IV. Irregular Code Optimization 

For a base code of high D stp but many degree 2 nodes, the 
lifted code does not generate satisfactory results since there are 
generally many sb with W sup = or 0.5, which benefit little 
from lifting. A degree augmentation (DA) method is provided 
to mitigate this mismatch between D s t p and M4up. 

• Base code optimization (DA+CA): Arbitrarily choose 
d u £ N. For those variable nodes of deg(x) > 3, 
construct d u (deg(x) — 2) pairs of auxiliary variable/check 
nodes (x a ,y a ) such that y a is connected to x and x a , 
Va € [1, d u (deg(x) — 2)]. Perform CA on the new graph. 
After CA, remove those auxiliary (x a ,y a ). 
m Perform cyclic lifting 

> Perform CA on the lifted code, focusing on the "local lift- 
ing sequence" instead of the "local edge configuration." 
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Fig. 3. Comparison among different CL constructions. "Typical CL" 
is randomly chosen from the CL ensemble using C4 as the base code. 
"CL+d2opt" is obtained by applying girth optimization over degree 2 nodes. 
By further applying CA, we have "CL+CA" which achieves D B t p = 12. 
"(DA+CA)+CL+CA" achieves D B t p = 13 and is obtained by optimizing 
the base code by (DA+CA). "Random Cnstr" (C5) and "Direct CA" (Ce) 
are plotted for comparison. The thicker and dashed lines are FER and BER 
curves respectively. 
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Fig. 4. Comparison between our codes in Fig. 3 and the best existing codes 
in [3], [5]. The thicker and dashed lines are FER and BER curves respectively. 



Let D' stp denote the stopping distance after DA+CA but 
before removing the auxiliary node pairs. For sufficiently 
large d u , the minimum W sup of the original graph grows 
linearly with respect to D' stp . The new base code thus has 
high Wsup and is suitable for cyclically lifting. The results 
of this DA-based method are plotted in Fig. 3. In Fig. 4, 
we compare the (DA+CA)+CL+CA code with some of the 
best existing randomized constructions in [3], [5], of which 
the asymptotic BEC thresholds are 0.07 and 0.22 worse than 
ours since their degree distributions are optimized for Gaussian 
channels instead. With similar error floors until 10~ 4 and 10~ 5 
respectively, our code has better threshold performance even 
after proper calibration. This discrepancy can be explained by 
the local vs. global edge rearrangement, and the latter of which 
has the potential of changing the threshold behavior. Both their 
degree distributions have 47-49% of variable nodes of degree 2 
while ours admits 62.5%. The large amount of degree 2 bits of 
our code leads to better threshold performance while limiting 



the lowest error floor one might achieve. In sum, our approach 
yields codes of low error floors comparable to the best existing 
codes without sacrificing the threshold performance even for 
codes of many degree 2 nodes. The CA applied in the last 
stage again eliminates any possible existence of performance 
outliers, a feature absent in the existing randomized construc- 
tions. Due to their distinct natures, numerical comparisons 
between the randomized CL construction and the algebraically 
constructed codes, e.g. quasi-cyclic codes [15], are deferred to 
a journal version of this paper. 

V. Further Research Directions 

Similar suppressing effect can be observed when consider- 
ing "trapping sets" for non-erasure channels. Exact formulas 
of the suppressing effect for the trapping sets and the cor- 
responding CA algorithm are currently under investigation, 
which utilizes the trapping set exhaustion algorithm in [16]. 

References 

[1] Y. Kou, S. Lin, and M. P. C. Fossorier, "Low-density parity-check codes 

based on finite geometries: A rediscovery and new results," IEEE Trans. 

Inform. Theory, vol. 47, no. 7, pp. 2711-2736, November 2001. 
[2] B. Vasic and O. Milenkovic, "Combinatorial constructions of low- 
density parity-check codes for iterative decoding," IEEE Trans. Inform. 

Theory, vol. 50, no. 6, pp. 1156-1176, June 2004. 
[3] X. Hu, E. Eleftheriou, and D. M. Arnold, "Regular and irregular 

progressive edge-growth Tanner graphs," IEEE Trans. Inform. Theory, 

vol. 51, no. 1, pp. 386-398, January 2005. 
[4] T. Tian, C. Jones, J. Villasenor, and R. Wesel, "Selective avoidance of 

cycles in irregular LDPC code construction," IEEE Trans. Commun., 

vol. 52, no. 8, pp. 1242-1247, August 2004. 
[5] A. Ramamoorthy and R. Wesel, "Construction of short block length 

irregular low-density parity-check codes," in Proc. ICC, 2004. 
[6] L. Dinoi, F. Sottile, and S. Benedetto, "Design of variable-rate irregular 

LDPC codes with low error floors," in Proc. IEEE Int'l. Conf. Commun., 

no. 1. Seoul, Korea, 2005, pp. 647-651. 
[7] E. Sharon and S. Litsyn, "A method for constructing LDPC codes with 

low error floor," in Proc. IEEE Int'l. Symp. Inform. Theory, July 2006, 

pp. 2569-2573. 

[8] C. Di, D. Proietti, E. Telatar, T. J. Richardson, and R. L. Urbanke, 
"Finite-length analysis of low-density parity-check codes on the binary 
erasure channel," IEEE Trans. Inform. Theory, vol. 48, no. 6, pp. 1570- 
1579, June 2002. 

[9] A. Amraoui, R. Urbanke, A. Montanari, and T. Richardson, "Further 
results on finite-length scaling for iteratively decoded LDPC ensembles," 
in Proc. IEEE Int'l. Symp. Inform. Theory. Chicago, 2004. 

[10] T. Richardson, M. A. Shokrollahi, and R. L. Urbanke, "Finite-length 
analysis of various low-density parity-check ensembles for the binary 
erasure channel," in Proc. IEEE Int'l. Symp. Inform. Theory. Lausanne, 
Switzerland, 2002, p. 1. 

[11] C. C. Wang, S. R. Kulkarni, and H. V. Poor, "Upper bounding the 
performance of arbitrary finite LDPC codes on binary erasure channels," 
in Proc. Int'l. Symp. Inform. Theory. Seattle, Washington, USA, July 
2006. 

[12] T. Richardson, "Error floors of LDPC codes," in Proc. 41st Annual 
Allerton Conf. on Comm., Contr., and Computing. Monticello, IL, 
USA, 2003. 

[13] S. H. Lee, W. H. Lee, J. J. Bae, S. I. Lee, and E. K. Joo, "Bit probability 
transition characteristics of LDPC code," in Proc. 10th International 
Conference on Telecommunications. Tahiti, French Polynesia, Feb. 
2003, pp. 553-557. 

[14] J. Gross, "Voltage graphs," Discrete Math., vol. 9, no. 3, pp. 239-246, 
Sept. 1974. 

[15] L. Chen, J. Xu, I. Djurdjevic, and S. Lin, "Near-Shannon-limit quasi- 
cyclic low-density parity-check codes," IEEE Trans. Commun., vol. 52, 
no. 7, pp. 1038-1042, July 2004. 

[16] C. C. Wang, S. R. Kulkarni, and H. V. Poor, "Exhausting error-prone 
patterns in LDPC codes," preprint. 



